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MIDCOURSE AND APPROACH GUIDANCE REQUIREMENTS 
FOR SIMPLIFIED ONBOARD CONTROL OF 
MOON-TO-EARTH TRAJECTORIES 

By Harold A. Hamer and Katherine G. Johnson 
Langley Research Center 

SUMMARY 

Simplified methods have been developed for midcourse and approach guidance of 
moon-to-earth trajectories. The methods utilize precalculated data to a great extent and 
require only simple onboard calculations based on optical angular measurements. These 
methods lead to guidance predictions sufficiently accurate for emergency control of entry 
angle. The approach guidance procedure can be adapted also to control the trajectory 
from earth-based line -of -sight measurements. 

An error analysis, for which a 10-second-of-arc measurement error and a 
0.2 -m/sec maneuvering error were assumed, has shown that an entry corridor of ±1° 
can be generally attained by use of approach guidance, with or without onboard midcourse 
guidance. Including the midcourse guidance correction insures higher accuracy at entry. 
The effects of measurement-star location and approximation errors on the guidance 
accuracy are discussed also. 


INTRODUCTION 

In Itmar missions navigation and guidance is normally accomplished by automatic 
procedures, which employ earth-based radar measurements. The inclusion of backup 
procedures which utilize onboard measurements is a desirable feature for emergency 
conditions which could occur during manned flights, such as failure in the radar system 
or loss of communications. 

In previous years various onboard guidance procedures have been developed for 
controlling moon-to-earth trajectories. (For example, see refs. 1 and 2.) Error 
analyses have shown these procedures to be sufficiently accurate; however, they generally 
require repeated optical measurements and a number of guidance maneuvers. The proce- 
dures presented herein differ in that the guidance correction is determined from a single 
onboard position fix made at a preselected time. When used in conjunction with certain 
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approximations derived from two-body theory, this position fix is adequate for either 
midcourse or approach guidance and uses a single guidance maneuver in each case. 

The guidance procedures presented herein represent a continuation of the work 
published in reference 3, wherein simplified guidance procedures were applied to earth- 
to-moon trajectories. For the most part, this paper incorporates the terminal guidance 
method of reference 3 to control entry angle of moon-to-earth trajectories. The method 
is employed as an approach guidance procedure either to refine the effects of transearth 
midcourse errors or to correct dispersions caused by transearth injection errors at 
the moon. In the latter case, no midcourse correction is made prior to the approach 
guidance correction. Similar measurements and guidance calculations are required for 
both cases, but there are differences in the accuracy characteristics and in the optimum 
measurement-star directions. A midcourse procedure which employs a fixed-time-of- 
arrival guidance law is also discussed, and some preliminary results are included. 

The accuracy characteristics of the methods were determined by use of a Monte 
Carlo procedure. Spherically distributed injection errors at the moon were assumed in 
both the position and the velocity of the spacecraft. At midcourse, the error distribution 
was based on work done with translunar onboard midcourse guidance (ref. 4.). The 
standard Jet Propulsion Laboratory n-body trajectory program (ref. 5) was used to 
generate all trajectories required for the analysis. 

SYMBOLS 


A,B,C 

D 




6D 

AD 

E = rg cos Vg 


matrices in midcourse-guidance equations (appendix A) 
position deviation from nominal trajectory 

position deviation in direction of star I, II, and III, respectively 
(appendix A) 

change in magnitude of D from T^ = 9.5 hours to T^ = 10 hours 
incremental value of position deviation 


K,Kj^,K2,Kg constants 

i,m,n direction cosine of line of sight to star with respect to X-, Y-, and Z-axis, 

respectively 
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earth radius 

range to earth center 

nominal entry range 

range to moon center 

incremental range to earth, r^^j - rj^ 

difference from nominal range to moon (table II) 

/ 2 2 2 ^^^^ 

position deviation from nominal trajectory, (Ax + Ay + Az y 

time to nominal atmospheric entry time 

time to nominal perilune time (table HI) 

time of first midcourse position fix (appendix B) 

time from transearth injection (perilune) 

time of first midcourse guidance maneuver (appendix B) 

difference from nominal entry time (table n) 

(2 2 

velocity deviation from nominal trajectory, (Ax + Ay + Az / 
spacecraft geocentric velocity 
spacecraft selenocentric velocity 

difference from nominal geocentric velocity (table II) 
difference from nominal selenocentric velocity (table II) 
guidance velocity correction 
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X,Y,Z 

x,y,z 

Ax,Ay,Az 

x,y,z 

Ax,Ay,Az 

Xr^yrj^r 

a 

y 

>'e 

Are 

6 

e 

X 

/i 


rectangular right-hand axis system in which X-axis is in direction of 
Aries, XY -plane is parallel to earth equatorial plane, and Z-axis is in 
direction of celestial north pole 


position coordinates in X,Y,Z system 


off-nominal position component in direction of X-, Y-, and Z-axis, 
respectively - for example. Ax = x^ - Xj^ 

the notation (Ax) represents the vector < Ay 


velocity coordinates in X,Y,Z system 


off-nominal velocity component in direction of X-, Y-, and Z-axis, 
respectively - for example. Ax = Xa - x^ 


position coordinates in rotating Cartesian axis system in which Xj.-axis 
lies along earth-moon line, X^Y^-plane is in earth-moon plane, and 
Zj.-axis is in northerly direction 


semisubtended angle 


flight-path angle at time of approach guidance 

flight-path angle at nominal entry altitude of 121.92 km (400 000 ft) 
incremental flight-path angle at nominal entry altitude, „ 

Q fO, 0^il 

angle formed at vehicle between line of sight to star and its projection in 
selenocentric orbital plane 


angle formed at vehicle by line of sight to star and line of sight to celestial- 
body center (referred to as star-to-body angle) 

angle between approach-guidance velocity vector and spacecraft velocity 
vector 


product of universal gravitational constant and mass of earth 
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cr 


standard deviation or root-mean-square value of error 



standard deviation or root-mean- square value of 


standard deviation or root-mean-square value of 


<^5 


submatrices in state -transition matrices 


s 

u 

(appendix B) 


Subscripts: 

a actual value 

D position deviation 

F first midcourse correction (appendix A) 

inj transearth injection 

m measured value 

meas due to measurement error 

n nominal value 

R earth radius 

r range to earth center 

S second midcourse correction (appendix A) 

T,pf at time of first midcourse position fix (appendixes A and B) 

T,1 at time of first midcourse guidance maneuver (appendix B) 

AV guidance velocity 

a semisubtended angle of earth 

y,e flight-path ar^le at nominal entry altitude 



e 


star -to -body angle 


1 to 6 specified directions 

A bar over a symbol indicates a vector. 

BASIC METHOD 

Two guidance techniques are considered in the present analysis: An approach 
guidance procedure is discussed fully, and a fixed-time -of -arrival midcourse guidance 
procedure is given in appendix A. Each method jdelds a guidance correction from a 
single position fix, which is determined from deviations in preselected directions from a 
nominal trajectory. For the midcourse guidance correction, three deviations in the 
direction of each of three stars are required, inasmuch as the spacecraft is guided to a 
fixed point. The approach guidance correction requires the determination of only one 
deviation. 

The geometry for the approach guidance method is shown in figure 1. The same 
variables are involved for the midcourse guidance. The deviation D of the actual 
trajectory from the nominal is normally determined by two onboard optical measurements: 
a star-to-body angle and a range measurement r^^j. Rai^e can be determined 

optically by measuring the angle subtended by the earth or by other measurements. (See 
ref. 4.) The corresponding nominal values would be known. 

The nominal trajectory chosen for this paper is depicted in figure 2. The perilune 
is considered as the transearth injection point and has a radius of 3358 km. The seleno- 
centric velocity is 1.992 km/sec; the time to earth atmospheric entry is 3.22 days. The 
nominal entry altitude selected was 121.92 km (400 000 ft). 

From an equation derived in reference 3 for control of periapsis radius, a guidance 
velocity correction is developed from a closed-form expression relating entry and 
upstream conditions. These conditions are then related to deviations from a nominal 
trajectory, and a technique is developed to control the entry angle y^. Although not 
specifically designed to control the along-track and cross-track position errors at entry, 
the method adequately compensates for these errors. 

For earth-approach trajectories, the substitution of the expression E = re cos 
in the aforementioned equation results in the following expression for correcting entry 
angle : 


6 



( 1 ) 



The alternate signs of the second term correspond to correcting for either side of the 
earth. The sign which results in the lesser value of AV would ordinarily be chosen. 

In equation (1), the flight -path angle y has the major influence on the magnitude 
of AV for a given guidance pointing angle X and guidance-maneuver time. It will be 
shown in subsequent figures that the value of y is highly correlated with the deviation 
from the nominal trajectory in certain directions. Furthermore, it will be shown that a 
measurement star in one of these preselected directions can be used to obtain the 
deviation D, and hence the required value of AV. 

APPROACH GUIDANCE PROCEDURE 

Application of the approach guidance procedure in two different ways was studied. 
First, the approach guidance correction is applied after a transearth midcourse correc- 
tion has been applied; second, the approach guidance correction is applied without prior 
midcourse correction. In essence, the first method corrects the approach trajectory for 
errors incurred at midcourse, the second corrects for errors incurred at the moon. 

With Midcourse Guidance 

It was assumed that a transearth first midcourse maneuver, employing the onboard 
method of reference 4, was performed about 20^ hours from perilune. (See fig. 2.) By 

lU 

using results obtained with this method for translunar trajectories, it was established that 
the midcourse guidance error is due mainly to errors in the onboard measurements, 
rather than in maneuver execution, and that the dominant measurement errors and the 
position-determination errors are generally in the direction of the spacecraft velocity 
vector. Thus, the midcourse guidance error is essentially an error in the magnitude 
of AV rather than in its direction. The resulting error distrubution yields an elongated 
ellipsoid with a la error of about 1 m/sec in the midcourse velocity correction. The 
la in-plane and out-of-plane errors in the direction of AV are 1.80° and 0.69°, respec- 
tively, where the reference plane is the nominal orbital plane of the spacecraft. For 
translunar trajectories, an error of 1 m/sec in midcourse AV resulted from a rai^e- 
measurement error of about 20 km. (See ref. 4.) From the equations presented in 
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appendix B, it was determined that for transearth trajectories, an error of 1 m/sec 
in AV was induced by an error of 40 to 50 km in determining range. 

By use of a Monte Carlo procedure, simulated midcourse errors were applied to the 
nominal trajectory to produce a number of perturbed trajectories. This method elimi- 
nates the calculation and application of both a first and a second midcourse correction. 

A second midcourse maneuver is inherent in fixed -time-of-arrival guidance, since it must 
correct the velocity error induced at the aim point by the derivation of the first midcourse 
velocity correction. Under actual operations, the approach-guidance measurements can 
be made only after a second midcourse correction has been applied. 

Star selection .- The deviation D was determined by use of the relationship 
D = Z(xa - xn) + m(y^ - y^) + n(z^ - 

where I, m, and n, the direction cosines of the line of sight to the measurement star,, 
are known. The measurement star is selected to yield the most effective direction for the 
deviation. Some of the directions for D which were investigated are illustrated in 
figure 3. Angles indicated represent attainable accuracies. The directions were chosen 
in two planes: the nominal orbital plane and the nominal instantaneous earth-moon- 
vehicle plane. At Tg = 9 hours, the angle between these planes is about 55°. Also, the 
spacecraft velocity vector V is about 13° out of the earth-moon-vehicle plane. As 
indicated, AV is generally preset 90° from the nominal velocity vector. The major 
axis of the position error ellipsoid, also indicated in the figure, lies within 1° of the 
orbital plane and was derived from errors incurred at first midcourse guidance. The 
deviation directions were chosen to essentially cover the entire spectrum of possible 
angles, as shown in figure 3(b). ^The deviation Dg is referenced to an earlier time 
Tg = 9.4 hours.^ 

The variation of deviation D with flight-path angle is shown in figures 4 to 7 for 
several directions of D and for several measurement times. Each data point repre- 
sents the condition at the indicated time on a perturbed trajectory resulting from errors 
at first midcourse guidance. Small differences in values of flight-path angle at the time 
of the measurement are magnified considerably at entry if not corrected: Large 
y values yield high g trajectories, whereas small values lead to trajectories which 
could skip out of or even miss the atmosphere. 

The data in figures 4 to 7 show that with midcourse guidance included, the devia- 
tion D is a good prediction of y. This relationship suggests that D can be used 
empirically to determine the guidance velocity correction because the AV required to 
change entry angle is dependent on the value of y. (See eq. (1).) It will be shown that 
this dependence is not as great for trajectories in which midcourse guidance is not 
included. 
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The scatter of the data in figures 4 to 7 produces error in the guidance procedure. 

It is important that the amount of scatter be minimized by selecting the most effective 
direction for D. The similarity of scatter characteristics in figures 4 to 7 suggests that 
for the case in which a midcourse guidance correction has been previously applied, the 
direction of the deviation (star) and the time of measurement are not critical so far as 
approximation error is concerned. This contrasts with limar-approach guidance where 
the direction of the star must be within 2° or 3° of a given plane. (See ref. 3.) For earth 
return, the trajectory position errors are concentrated along the major axis of the error 
ellipsoid (fig. 3); for example, at Tg = 9 hours, the length (la) of the major axis is 
236 km compared with 9.4 km and 2 km for the other axes. Similarly, the velocity errors 
are almost entirely along the major axis of the velocity error ellipsoid, which is in the 
approximate direction of the earth. 

It should be pointed out, however, that the deviation direction affects the measure- 
ment sensitivity which in turn affects the guidance accuracy attributable to measurement 
error. Comparison may be made between figures 4(b) and 7(a) where the deviations Dj 
and D 3 are both perpendicular to the range vector but in a different plane. Choosing 
the deviation in the earth-moon-vehicle plane (fig. 7(a)) reduces the deviation (measure- 
ment) sensitivity considerably, and the reduced sensitivity increases the effect of 
measurement error. 

The results in figures 5 and 6 show deviation directions which yield maximum 
sensitivity. Figure 6 gives the variation of flight-path angle with Ar, the difference 
between the nominal and measured range. The Ar values essentially indicate devia- 
tions measured directly toward (or away from) the center of the earth. 

The deviation directions in figure 5 would lead to large measurement-error effects 
on guidance accuracy since 9 ^ 90°, and even larger effects in figure 6 where = 0°. 

To minimize measurement error, the most suitable deviation direction is perpendicular 
to the nominal range vector, as shown in figures 4 and 7(a). These results are discussed 
in the section "Approach -Guidance Accuracy Characteristics." Figures 7(b) and 7(c) are 
included to give broader coverage of the entire spectrum of directions considered in the 
earth-moon-vehicle plane. (See also fig. 3(b).) 

The foregoing results for approach guidance have been applied to the case in which 
midcourse guidance velocity errors are generally along the direction of the spacecraft 
velocity vector, the direction normally expected for onboard midcourse velocity errors. 

It is of interest to examine the effect on approach guidance of midcourse errors in other 
directions which might result from other midcourse procedures. Accordingly, midcourse 
velocity errors of the same magnitude and angular displacement previously used were 
applied normal to the spacecraft velocity vector in the orbital plane. The resulting data 
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are presented in figure 8. (Note the staggered vertical scale, which is read in such a 
manner that at D = 0, y = -73.954°.) Although a strong correlation between flight -path 
angle and certain directions of deviation D is shown, the range of effective directions 
was restricted. For example, the scatter for D5 is considered imacceptable in that the 
maximum error in y is >0.01°. The deviations D2, D3, and D4, and Ar also 
exhibited excessive scatter, as well as low sensitivity, and are not shown. It is apparent 
that should the midcourse velocity errors be perpendicular to the spacecraft velocity 
vector, satisfactory deviations can be tak'in only in a limited range of directions. (See 
fig. 3(a).) 

Guidance velocit^requirements.- Characteristics of the approach guidance velocity 
are shown in figures 9 to 13. Figure 9 is an example of the variation of the velocity 
correction with the deviation D3. The positive and negative values of AV indicate 
A = ±90°. The circular symbols represent various trajectories perturbed at first mid- 
course guidance and were computed from the curve presented in figure 7(a) for 
Tg = 9 hours. Some data are also shown for two other times to indicate the general 
trend of the data. The offset in AV at D3 = 0 (for each curve) is the results of two- 
body approximation (eq. (1)) and must be corrected by shifting the curves vertically to 
zero offset. After obtaining the value of the deviation, the astronaut uses a curve such as 
that shown in figure 9 to determine the required AV magnitude. 

The calculations for AV are performed before flight and, as shown by equation (1), 
require the perturbed values of y, r, and V at the measurement time, as well as the 
nominal value of E = rg cos y^. The AV values in figure 9 represent application of the 
thrust in the nominal orbital plane, perpendicular to the nominal velocity vector, the 
optimum direction for most times along the approach trajectory. (See fig. 10.) 

Given in figure 10 are the magnitudes of AV required for correcting one particular 
perturbed trajectory; any other trajectory would have similar requirements percentage- 
wise. The astronaut uses this type of data to choose the time for the guidance maneuver. 
The data show that the AV requirement increases as the guidance maneuver is delayed 
to times closer to the earth. The value of the guidance pointing angle A = 90° represents 
application of the AV vector perpendicular to the nominal velocity vector; the value 
A = 0° represents application of the AV vector along the nominal velocity vector. The 
AV vector is always applied in the nominal orbital plane. It is evident that A = 90° is 
essentially optimum for all times from entry. For T^ = 4, the optimum value is about 
70°, but the use of A = 90° results in a negligible increase in the AV requirements. 

The results in figures 11 to 13 define the approach guidance velocity requirements. 
The symbols in figures 11 and 12, which correspond to the data of figure 9, have been 
shifted for zero offset. Note in figure 11 that only three of the 50 perturbed trajectories 
exceed the Ay^ value of ±1°, and these by only a small amount. Thus with midcourse 
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guidance included, approach guidance becomes less important except when the midcourse 
error is not along the spacecraft velocity vector. In figure 8 , where the error is essen- 
tially normal to the velocity vector, the magnitude of Ay^ exceeds ± 1 ° two-thirds of the 
time. As will be discussed subsequently, data in the form shown in figure 11 are used to 
determine the effect of scatter (approximation) error and guidance velocity -cutoff error 
on entry-angle accuracy. Equivalent data are shown in figure 12. The AV values are 
shown as a function of rg cos because this relationship can be calculated analytically 
from the derivative of equation (1). (See appendix C.) The calculated results and those 
obtained from the slopes of the curves in figure 12 coincide with the curve shown in 
figure 13. The curve shows the substantial increase in the guidance velocity requirement 
as the guidance maneuver is delayed to times closer to the earth. Obviously, the velocity- 
requirement ratio A(AV)/A^re cos y^ can be readily converted to the more usable 
form AV^Ay^ for error analysis. 

Sensitivity o f entry angle to guidance measurement .- Examples of the variation of 
deviation with incremental entry angle are shown in figure 14. This information is 
required for determining the effect of measurement error on the control of entry angle. 
The data are shown with respect to the deviation rather than D 3 since the effect 

of measurement error for the Dj direction is smaller because of the increased sensi- 
tivity of D. 


Without Midcourse Guidance 

The approach guidance method discussed in this section corrects perturbed trajec- 
tories due to trajectory errors incurred at the moon and is the only guidance applied to 
the return trajectory. For the analysis, the Monte Carlo procedure was used to produce 
a number of perturbed trajectories emanating from perilune. The errors were essen- 
tially distributed spherically with la values of about 2 km for position deviation s and 
about 2 m/sec for velocity deviation u. These errors could represent errors in the 
transearth injection for return to the earth. 

Star selection .- Shown in figure 15 are examples of the variation of y with D 
at Tg = 9.4 hours for two star directions. (Note the staggered scale, which is read in 
such a manner that at D = 0, y = -74.11°.) Except for the ranges covered, the data are 
similar to those shown in figures 4 to 7, where midcourse guidance was included. It is 
apparent that under certain conditions, the deviation D essentially predicts y. 
However, as shown by figures 15 and 16, the deviation direction which gives the best cor- 
relation for y is not the best for predicting the required approach correction AV. 

The variation of AV with deviation is shown in figure 16 for the perturbed trajec- 
tories. This figure is similar to figure 9, except that certain perturbed trajectories do 
not enter the atmosphere if no AV correction is applied. The AV values in figure 16 
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were calculated from the actual y values of figure 15. The rms value of AV shown 
for each set of data in figure 16 is the Ict value of the scatter in the data about the line. 
This value is used in the subsequent error analysis as a measure of the amount of scatter 
(approximation) error. As indicated in figure 16, the measurement- star direction is of 
great importance in case of no midcourse guidance. In regard to scatter, the optimum 
direction of D at Tg = 9.4 hours corresponds to a star lying in the nominal orbital 
plane in a direction 77° from fn- (Note in fig. 3(a) that this direction is to the left of the 
earth.) 

The variation of AV with D for the perturbed trajectories is shown in figure 17, 
for two values of 6^ ^.t Tg = 17.4 hours. Here again, the Ic values correspond to the 
scatter about the lines. (Note the staggered vertical scale, which is read in such a 
manner that at D = 0, AV = -0.6 m/sec. The curves should be shifted vertically to 
correct for this zero offset.) For Tg = 17.4 hours, it is apparent from the scatter that 
the optimum direction of D corresponds to = ^0°. Data are shown for 6^ = 90° 
because no range measurement is required. Elimination of the rai^e measurement is 
important from an operational standpoint, as well as for the fact that range-determination 
error is a dominant factor affecting the guidance measurement accuracy. The effect of 
the increase in scatter incurred by omitting the range measurement is discussed in the 
subsequent error analysis. 

With regard to scatter error, the optimum direction of D at Tg = 4.4 hours is 
about 71°. It is of interest to note that since the true anomaly of the trajectory changes 
approximately 15° from Tg = 17.4 hours to Tg = 4.4 hours, the corresponding change in 
the optimum direction of the measurement star is 6° away from the earth. It should be 
noted also that for acceptable scatter, the nominal value selected for the star-to-body 
angle may lie within a region of several degrees in either the in-plane or out-of -plane 
direction. 

Guidance veloci ty req uireme nts.- The gmdance velocity requirements, with respect 
to the deviation D, were discussed in the previous section in connection with selecting the 
optimum measurement star. In figure 18, AV is shown as a function of Ay^ for 
T„ = 9.4 hours, where the quantity Ay is the difference between the nominal and actual 
values of entry angle (at an altitude of 121.92 (400 000 feet)) for the various perturbed 
trajectories. The AV values are those given in figure 16 and are calculated values 
required to correct the Ay^ to zero. There are perturbed trajectories beyond the point 
Ay = 4°, but these do not enter the atmosphere. The sensitivity of AV with Ay 

G G 

at Ay^ = 0 is utilized in the error analysis. 

Sensitivity o f entry angle to guidance measurement.- The values in figures 16(b) 
and 18 were cross plotted to obtain figure 19, which shows the variation of D0 with Ay^. 
The value of the slope of this curve at Ay^ = 0 is essential in determining the effect of 
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measurement error on guidance accuracy. Only the slope at Ay^ = 0, that is at nom- 
inal y^, is required because all perturbed trajectories are corrected to the vicinity 
of Ay^ = 0. It is of interest to note that the slope in figure 19 is about twice that for 
data where the trajectories were perturbed at first midcourse guidance (fig. 14). Absence 
of scatter in the data shown in figure 19 can also be used as a good indication for the opti- 
mum direction of D. 

APPROACH- GUIDANCE ACCURACY CHARACTERISTICS 

In this section, the more important errors associated with the approach guidance 
procedure are defined and analyzed. The analysis spans the region for performir^ the 
guidance from Te = 18 hours to Tg = 4 hours. At Tg = 18 hours the spacecraft is 
about midway between the earth and moon. (See fig. 2.) Both procedures, that is, with 
and without midcourse guidance, are analyzed. 

The approach guidance procedure is designed to control only the entry flight-path 
angle; however, the error analysis has shown that the along-track and cross-track posi- 
tions are also controlled to a reasonable accuracy. For example, in the case with a 
midcourse guidance correction included, the approach guidance procedure reduced the 
la position error at entry from 122 km to 17 km, these errors being almost entirely in 
the along-track direction. Similar results were obtained for the case without a midcourse 
guidance correction, except the uncorrected position errors at the earth were considerably 
higher. 


Effect of Measurement Error 


Error in the onboard optical measurements affects the approach-guidance accuracy. 
Error equations for such measurements were developed in reference 6. The equation 
used in this report is 


OD = 



(r sin B a q) 



( 2 ) 


which corresponds to uncorrelated errors in the measurement of range r and star-to- 
body angle B. It is assumed that range is measured by a, the half-angle subtended by 
the earth at the spacecraft. The error Oq is constant at all values of time; whereas Oj. 
varies with time (or range) according to the relation given in reference 6, 

1/2 

(3) 


°^r = 
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where ctj^ is the uncertainty in defining the earth's radius and <j^ is constant. The 
variations of the nominal range and a with time are shown in figure 20. 

In equation (2) nominal values for r and 6 are used. For values of 6^ « 90°, 
it is seen that 


In this case, the range-measurement error is insignificant; hence, the nominal range 
value can replace the measured range in the equation given for D in figure 1. Because 
no range measurement is required, the approach method could be adapted to control the 
trajectory from earth-based line-of-sight measurements. If, for instance, a failure in 
the radar system prevented range and range-rate measurements, earth-based angular 
measurements could be substituted for star-to-body measurements. Even though rela- 
tively inaccurate, the angular measurements could be averaged over a period of time. 

For ^ 90°, the effect of range-measurement error Or becomes important, as 
shown by equation (2). Figure 21 shows the variation with time of Oj., the la error in 
determining range from the a measurement. The data are shown both with and with- 
out consideration of the earth-radius uncertainty aj^. The effect of excluding this 
uncertainty in the ensuing error analysis would be negligible. As shown by the other 
curve in the figure, it may be more accurate to use the nominal value instead of measuring 
the range at certain distances from the earth. This condition would, of course, depend 
upon the injection-error statistics at the moon. The curve was determined from normal- 
ized values obtained from the perturbed trajectories. The value of the injection -position 
error Og is given in the figure although its effect on a^. is minor. 

The effect of measurement error on the guidance accuracy is presented in figure 22 
for approach guidance with and without a previous midcourse correction. The upper plot 
shows the variation of error in determining the deviation D with time from entry. 

These data were obtained by the use of equation (2) and the solid curve in figure 21. The 
sensitivity ratios shown in the middle plot of figure 22 are used to convert the Oq data 
to the rms values shown in the lowest plot. Sensitivities were obtained from data 
such as those shown in figvures 14 and 19. For the data with midcourse guidance, the sen- 
sitivity ratio is nearly constant with time since the overall magnitude of D does not 
chaise appreciably. (See figs. 16 and 17.) 

The effect of errors in the range measurement is shown in the two curves in the 
upper plot of figure 22. The upper curve was calculated for the values of which 
correspond to minimum scatter. (For example, see figs. 16(b) and 17.) Since star direc- 
tion is not critical for scatter in the case with midcourse guidance included, the lower 
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curve was calculated for 0^ = 90°, which is optimum with regard to measurement error. 
The deviation error is highly dependent on the value selected for 9 ^. For instance, 
for D 2 , which corresponds to » 126.3°, Op at Tg = 9 hours would be about 55 km. 
K Ar were used as the guidance measurement, 0^ ~ 9°, (fig- 6), the value of ctq would 
equal that of ay, which is about 90 km at that time. In these cases, however, the values 
of Cy g would not increase in the same proportion because the sensitivity ratio Ay^y^AD 
would be about three times smaller in the case with midcourse guidance. 


Effect of Approximation Error 

The upper plot in figure 23 shows the effect of scatter (approximation) error on the 
accuracy of controlling entry ai^le. These curves were determined from 



scatter 



scatter 



(4) 


where the value for was obtained from data such as those shown in fig- 

^ / scatter 

ures 16 and 17. The ratio Ay^yAV is the reciprocal of the velocity-requirement ratio, 
which can be determined from figure 13 by use of the equation 


AV ^ a(av) 
Ayg A(rg cos yg) 


^^^(0 051 ^^ ( 0 . 0001 ) 


where 65 083 000 meters is the nominal entry range and 0.051° is the change in yg from 
its nominal value of -6.281° due to a change of 0.0001 in the cosine function. 

The value of (*^y,^ggg^j-j.gj. was also checked at several times aloi^ the trajectory 
by simulated guidance corrections. For example, the perturbed trajectories in figure 9 
were corrected by usii^ the faired values of AV at Tg = 9 hours. Each trajectory was 
then propagated to the nominal entry altitude. The rms value of the entry -angle error was 
determined and was found to compare closely with that calculated by equation (4). 


Effect of Maneuvering Error 

The velocity -requirement ratio is also employed to determine the effect of guidance 
maneuvering error. Velocity -cutoff error comprises the major portion of this error. 

(See ref. 6.) 

The rms entry-angle error due to velocity-cutoff error is shown in figure 23 for a 
typical IcT value of 0.2 m/sec. This cuiwe was determined by multiplying the reciprocal 
of the velocity-requirement ratio by 0.2 and corresponds to the case in which the guidance 
pointing angle X is 90°, which is optimum for AV magnitude. It should be pointed out 
that by sacrificii^ the requirement for the AV magnitude, which is small — especially 
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in the case where midcourse guidance has been previously applied - and applying the 
AV vector at X = 0°, the effect of the maneuvering error can be reduced by about two 
thirds. This same effect could be achieved by delaying the guidance maneuver. To 
obtain the two-thirds reduction, however, would require considerable delay (fig. 10), 
which may not be operationally feasible. 


Combined Effect of Guidance Errors 


The lowest plot in figure 23 shows the combined effect on entry-angle error due to 
measurement error, scatter error, and velocity-cutoff error according to the equation 

-il/2 


y,e 




\2 


+ (O'. 


\2 


vmeas ^ scatter 




cutoff 


With midcourse guidance .- Statistically speaking, the effect of measurement error 
has only a minor role in the overall guidance accuracy when a midcourse correction has 
been included; the velocity -cutoff error has the major effect. These facts are apparent 
when the curve for total a is compared with the curves for these two effects shown 

7,6 

in figures 22 and 23. It is significant to note that even if the instrument measurement 

error were doubled to 20 seconds of arc, the overall a would be little affected. If 

■y,e 

the velocity -cutoff error were reduced in the manner previously described, the total 

a would be negligible. 

7,e 

The data shown in figures 22 and 23 pertain to midcourse velocity errors, generally 
in the direction of the spacecraft velocity vector. An error analysis was performed for 
the situation where these errors were approximately normal to the velocity vector (fig. 8), 
and the results were found comparable with those shown for total a in figure 23, 

n . y 9^ 

providing = 90 . The main difference is that Ay^/AD more than doubles, and there- 
fore the effect of measurement error increases. 


Without midcourse guidance .- For the method without midcourse guidance compar- 
ison of the curve showing total effects and curves showing the three separate effects indi- 
cates that error in entry angle for this method is due chiefly to measurement error. 

If 6 = 90° at Tp = 17.4 hours, the total a would be due mostly to scatter (fig. 17) 
n y 

and would be about the same as that shown in figure 23 where = 80^. If the curve in 
figure 21 representing no range measurement (r^^ = r^j) could be applied, the error in 
entry angle could be reduced for times earlier than Tg = 12 hours. 

Maneuver time .- In regard to the total accuracy characteristics shown in figure 23 
and for an entry corridor of ±1°, it is apparent that with midcourse guidance, the approach 
guidance can be selected at any time along the trajectory because the 3cr value is always 
below 1°. For approach guidance without midcourse guidance, it appears from the errors 
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shown that an approach maneuver at about Tg = 8 hours or later would result in 
3a values of error below 1°. At Tg = 8 hours the average AV requirement is 
about 5 m/sec for the method without midcourse guidance and only about 0.5 m/sec for 
the method with midcourse guidance. The latter procedure, however, requires an addi- 
tional average AV of about 3 m/sec for the midcourse maneuver under normal 
conditions. 


CONCLUDING REMARKS 

Results were obtained for simple onboard guidance of moon-to-earth trajectories 
which can be applied either with or without a midcourse correction. Normally one, but at 
most two, onboard angular measurements are required, from which a simple calculation 
yields the magnitude of the guidance velocity required to correct the trajectory. It was 
shown that the method which included the onboard midcourse correction is far superior. 
This method has much higher accuracy in controlling the entry angle and has more flex- 
ibility in the star selection. The method without midcourse guidance could be resorted to 
under the extreme condition that a midcourse correction were not available on the return 
trajectory. 

Under certain conditions, a rar^e measurement is not required; hence, the method 
could be applied to a procedure whereby the approach trajectory is controlled by earth- 
based line-of-sight measurements. 

An error analysis showed that errors in the onboard guidance measurements do not 
primarily affect the overall approach-guidance accuracy if midcourse guidance has been 
included. Furthermore, the effect of guidance maneuvering error can be reduced by as 
much as two-thirds by changing the direction of the thrust vector or by delaying the 
guidance maneuver. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 8, 1971. 
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APPENDIX A 

SOME NOTES ON FIXED-TIME-OF -ARRIVAL ONBOARD MIDCOURSE 
GUIDANCE FOR TRANSEARTH TRAJECTORIES 


One objective of the approach guidance procedure given in the main text of the 
report is to refine the effects of midcourse guidance error. Hence, the use of some type 
of midcourse guidance procedure was required as a prerequisite for the analysis of the 
approach procedure. The onboard method of reference 4 was employed for this analysis. 
In this reference, the procedure for calculating the midcourse maneuver velocity correc- 
tion was developed from standard guidance equations, which make use of transition-matrix 
theory. The well-known fixed-time-of-arrival guidance, which guides the spacecraft to a 
given aim point on the nominal trajectory, was employed. Two basic equations which 
concern the navigator are 


I 


j 


AVp = [a] <i 



1 *^2% > 



T,pf * 





AVg . [C] 


K,D, 
^2°II ) 




T,pf 


The second midcourse correction AVg corrects the velocity error induced at the aim 
point by the derivation of the first midcourse velocity. The quantities A, B, and C 
are 3-by-3 matrices which are precomputed from nominal state -transition matrices and 
the direction cosines of three stars used for the measurements. The deviations D are 
the three position components as measured by the deviations of the perturbed trajectory 
from the nominal trajectory in the directions of the three stars. The constants K are 
predetermined from the variations shown in figure 24. The relation is 


The values in figure 24 were obtained from a sample of trajectories randomly perturbed 
at transearth injection. (See table I.) Figure 25 shows the geometry for the angles and 
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APPENDIX A — Continued 


directions referred to in figure 24. The time = 9.5 hours corresponds to the time of 
the position fix, and the time Tj. = 10 hours is the time of the midcourse guidance 
maneuver. 


Accuracy Characteristics 

In reference 4 the method was thoroughly studied and found to be sufficiently 
accurate for midcourse control of translunar trajectories. Some results are presented 
herein which indicate that the method would apply as well, or better, to transearth 
trajectories. 

Effe ct of star selection .- One important characteristic of the method is the selection 
of the three measurement stars in directions such that 6D can be accurately predicted 
from D calculated for the time of the position fix. This procedure eliminates the need 
for a second position fix and, hence, eliminates much of the effect of measurement error. 
For translunar trajectories, the region of acceptable stars for measurement is about 40°, 
as shown in sketch (a). 



Suitable 

stars 


Figure 24 shows that for transearth trajectories, good prediction of 6D may be obtained 
for an even larger region of stars. The only unacceptable direction is the bottom curve, 
characterized by large scatter and small sensitivity in D. Figure 25 and the upper five 
curves in figure 24 show that the acceptable region covers an angle greater than 90° in 
the selenocentric orbital plane as well as large out-of -plane angles. The region diago- 
nally opposite this region would also be acceptable as may be seen in the sketch of the 
trajectory geometry (sketch (b)). (Values of 5 near 90° provide no information on 
change in the orbital plane and should be avoided.) 
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APPENDIX A - Continued 


To earth 



Sketch (b) 

Effect of trans ition -matrix theory.- The midcourse guidance method employs 
transition-matrix theory; therefore, another important characteristic of the method is the 
effect caused by the linear approximation made in using this theory. The state-transition 
matrix is strictly applicable only when the equations of vehicle motion are linear. 

Because of the actual noiHinearity of these equations, the midcourse guidance method is 
limited to perturbed trajectories that are reasonably close to the nominal trajectory. 

The data in table II show the results of correcting a sample of 33 perturbed transearth 
trajectories with the onboard midcourse method. No errors due to measurement error 
or in executing the guidance maneuver are included; the errors at the earth are due only 
to the effects of the linear approximation made in using the transition -matrix theory. As 
indicated by the large errors at the moon (at nominal perilune time), which could corre- 
spond to the transearth injection errors, the perturbed trajectories are widely dispersed 
about the nominal. The large dispersions are even more apparent at first midcourse 
guidance, which is about 20.5 hours beyond perilune time. Such large perturbations for 
injection out of lunar orbit are extremely improbable; however, these perturbations could 
possibly apply imder emergency conditions when a translunar free-return trajectory is 
altered on approach to the moon to insure a 100-percent probability of miss. 

The last two columns in table II indicate that widely dispersed trajectories can be 
reasonably controlled by guidance employing transition-matrix theory. It is interesting 
to note that the errors at the nominal entry altitude are smaller than those at the aim 
point. It is also seen that without the addition of the second midcourse velocity correc- 
tion, the errors in spacecraft velocity and flight-path angle are small at the nominal 
entry altitude. These errors would not be reduced by including the second midcourse 
correction. 


20 



APPENDIX A - Concluded 
Velocity Requirements 

A comparison of midcourse-guidance velocity requirements of the onboard method 
for translunar and transearth trajectories is shown in table m. To establish the 
perturbed trajectories, injection errors of the same magnitude were applied at the earth 
and at the moon. These errors were applied by randomly changing components of the 
position by 10 km and components of the velocity by 10 m/sec. (See table I.) The velocity 
requirements pertain to a guidance maneuver time 10 hours from translunar or transearth 
injection. It can be noted in table III that the translunar requirements are 3 to 4 times 
higher than the transearth requirements. As expected, in both cases, the first midcourse 
velocity requirement is fairly insensitive to time selected for the aim point. However, 
the second midcourse velocity requirement does depend on the aim-point time, with the 
translunar (transearth) values decreasing (increasing) as the aim point is moved closer 
to the target body. The large transearth value for the aim point at nominal entry time is 
attributed to the relatively large magnitude of the spacecraft velocity at this point. 
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APPENDIX B 


ESTIMATION OF ONBOARD MIDCOURSE VELOCITY ERROR 

The magnitude of the guidance velocity error for onboard midcourse procedures can 
be estimated to a good degree of accuracy. The general form of the equation for onboard 
midcourse guidance, developed from reference 4, is 



Since | | is essentially a vinit matrix, then the following approximation can be made: 

(Bl) 

As shown in reference 4, error in the range measurement has the dominant effect on the 
error in AV. Also shown is the fact that for trajectories perturbed because of injection 
error, the position deviations Ax^ pj and Axrp are essentially in the direction of the 
range vector. These conditions permit the AV magnitude error to be determined by 
inserting appropriate errors for the position deviations in equation (Bl). To accomplish 
this, a standard-deviation value of the range -measurement error can be determined for 
the time Tpj, based on instrument inaccuracies, and then substituted for Axrp pj as 
the error vector (in the range direction). The error vector Axrp j, which is somewhat 
larger in magnitude, can be determined from any perturbed trajectory by the percentage 
change in Ar from the times Tp^ to T^. 
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APPENDIX C 


VELOCITY REQUIREMENT FOR APPROACH GUIDANCE 

An analytical expression for the variation of approach guidance velocity AV with 
the entry condition rg cos is presented in this appendix. The velocity-requirement 
ratio A(Av)y^(rg cos y^) is useful for error analysis. The equation which follows is 
the exact expression for this ratio (at small values of AV) and was determined by 
differentiating equation (1) with respect to rg cos y^. Except for distances very close to 
the earth, only the middle term of the equation is significant. For example, at 
Tg = 9 hours (r “ 110 000 km) the first and last terms contribute only about 1 percent to 
the total value. 


9(AV) 2r2vE cos (y + X) [cos y + cos X cos (y + X)[| 


[r^ cos^ (y + X) + gjn2 x + |r2 cos^ (y + X) - 


[t2 cos 2 {y + X) - E^ 


E2(v2sin2xy^-^ 
Pfs r 


[r^ cos^ (y + X) - e2] - e2) v 2 sin2 X + [r2 gos2 (y + _ 


where 


E = rg cos y^ 
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TABLE I.- TRANSLUNAR AND TRANSEARTH PERTURBATIONS USED 
TO DETERMINE MIDCOURSE GUIDANCE CHARACTERISTICS 
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TABLE II.- THE la RESULTS FOR WIDELY DISPERSED 
TRANSEARTH TRAJECTORIES BEFORE AND AFTER 
"PERFECT" MIDCOURSE GUIDANCE CORRECTION 


Quantity 

Before correction 

After correction 

At nominal 
perilune time 
(transearth 
injection) 

At first 
midcourse 
guidance 

At aim point 
(nominal 
entry time) 

At nominal 
entry altitude 

1 s, km 

896 

3323 

284 

65.3 

u, m/sec 

208 

41.8 

^243 

1 

Ar^, km 

125 

1558 


1 

Ar, km 



36.4 


6V^, m/sec 

27.7 

17.1 



6V, m/sec 



^39.8 

bl.33 

Ayg, deg 




bO.29 

At, sec 




28.8 


^ AVs included. 
^AVs not included. 
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TABLE in.- VELOCITY REQUIREMENTS FOR 


FIXED-TIME-OF-ARRIVAL GUIDANCE 



at 10 hours from 


inj ectio^ 


Location of aim point 

Average values 

AVp, m/sec 

iVs/iVp 

Translvinar trajectories 

Lunar sphere of influence ^Tp » 14.6 hr^ 

55 

0.13 

Nominal perilune time 

54 

0.08 

Transearth trajectories 

Tg = 17 hr 

15.7 

0.17 

Te = 9 hr 

15.2 

0.16 

Tg = 1 hr 

14.7 

0.27 

Nominal entry time 

14.7 

3.98 






(a) Nominal orbital plane. 


To moon 



(b) Nominal instantaneous earth -moon -vehicle plane. 


Figime 3.- Illustration showing direction of deviations studied for Tg = 9 hours. 
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(a) Deviation perpendicular to nominal geocentric range vector. 

Figure 7.- Flight-path angle correlation with deviations taken in instantaneous 
earth-moon-vehicle plane at Tg = 9 hours. 







Flight-path angle. 



Figure 8.- Flight -path angle correlation with deviations taken at Tg s 9 hours. 
Midcourse guidance velocity errors are approximately normal to spacecraft 
velocity vector. (Note staggered scale.) 
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! Figure 9.- Correlation of approach guidance velocity with deviation D3. 
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(a) = 80°. 

Figure 16.- Correlation of approach guidance velocity with deviation at Tg = 9.4 hours. 

X = ±90°: D is in nominal orbital plane. 



Guidance velocity, AV, m/sec 



D, km 

Figure 17.- Correlation of approach guidance velocity with deviation at Tg = 17.4 hours. 
X = ±90°; D is in nominal orbital plane. (Note staggered vertical scale.) 



Guidance velocity, AV, m/sec 



Figure 18.- Approach guidance velocity requirement as a function of incremental 
entry angle at nominal entry altitude. X = ±90°; Tg = 9.4 hours. 




Range -determination 



Figure 21.- Range-determination error. Measurement error = 10 seconds of arc. 
Dashed line includes effect of 2-km uncertainty in earth radius. 
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Time to nominal entry time, T^, hr 


Figure 23.- Approach-guidance accuracy characteristics with effects 
of scatter error and maneuvering error included. The la value of 
velocity-cutoff error is 0.2 m/sec. 




Figure 24.- Change in deviation from Tt = 9.5 hours to = 10 hours for various 
perturbed transearth trajectories. D and 6D agree in sign except for points 
below the zero line. 



To earth 



Figure 25.- Geometry in vehicle selenocentric orbital plane 
at 9.5 hours from transearth injection. Positive 6 is 
in northerly direction from plane; 6 is not necessarily 
in plane. 
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